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Interpolation Properties of Orthogonal Sets of Solutions 

of Differential Equations* 

By 0. D. Kellogg. 



As the present paper is of the nature of a continuation of the one pub- 
lished in the last number of this Journal, an extended introduction may be 
dispensed with. Suffice it to say that we shall here be concerned with the 
problem of extending to the orthogonal function sets arising from ordinary 
differential equations of second order, the properties there derived for sets 
arising from integral equations. 

1. The Differential Equation and Boundary Conditions. 
We shall be concerned with differential equations of the form 

;£(*<*) ^)+(fl(*)+M*))*=0, 

v ' p>0 for 0<as<l, r>0 for 0<*<1, (1) 

with the general homogeneous self -ad joint boundary conditions 

ay(l)+by'(l)=cy(0)+dy'(0), a'y(l) +b'y' (1) =Cy(0) +d'y' (0). (2) 

If p(x) vanishes at one of the end points of the interval (0, 1), but in 
such a way that the differential equation has a solution which remains finite 
and different from zero, with a finite derivative in the neighborhood of this 
point, one of the boundary conditions (2) is to be replaced by the condition: 
y (x) remains finite and different from in the neighborhood of this point. If 
p(x) vanishes at both end points, both equations (2) are to be replaced by this 
demand for both end points. Examples of these cases are: (a) the functions 
J (a x) y/oc, J^{a x x) V#, . . . ., where a Q , a ly . . . . are the successive roots of the 
Bessel function of zero order J (%); (b) the Legendre polynomials on the 
interval ( — 1, 1). 

* Presented to the American Mathematical Society, December 1, 1917. 
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226 Kellogg: Interpolation Properties of Orthogonal 

The existence of solutions of the problem (1) and (2), and their oscilla- 
tion properties have been extensively studied,* and we shall make use of some 
of the more common ones. Our interest will center rather in the interpolation 
properties, and some of the consequences of the property (2>),t which it will 
be our object to establish. 

By "harmonics" we shall understand solutions of the differential equation 
and boundary conditions. The corresponding " frequencies " are the values of 
Jl for which these solutions are possible. There is an infinite sequence of 
frequencies, without finite limit point, and bounded below. We shall think 
of them written in ascending order of magnitude, a frequency being written 
twice if it corresponds to two harmonics. The harmonics are orthogonal ; t 
even two corresponding to the same frequency may be made orthogonal by a 
proper choice of the integration constants entering them. No harmonic van- 
ishes simultaneously with its derivative. The zeros of any harmonic separate 
those of any harmonic with the same or lower frequency. In addition to these 
known facts, we assume explicitly (1) the continuity of the coefficients of the 
differential equation and of the derivatives involved, (2) that the boundary 
conditions are such that for any pair of harmonics 



P(0) 



fc(0),$,(0) 
*I(0),*J(0) 



--p(l) 



fc(i),<Mi) 



(3) 



and (3), that the i-th harmonic has exactly i zeros in the interior of (0, 1) for 
all i. A harmonic is " even " or " odd " according as the number of interior 
zeros is even or odd. 

2. Types of Boundary Conditions. 

For what follows it will be convenient to separate the boundary conditions 
into types according to the behavior of the corresponding harmonics at the 
end points of (0, 1). Consider first the case in which the determinant of 
coefficients in (2) ab' — a'b vanishes. There will then be a relation of the form 
a"<M0) +&"<?>«(0) =0. From this it follows that the function 

M(x)=p(x)[<p' i (v)<l> j (x)-cl> i (x)ti(x)] 

* See B6cher, '< Legons sur les methods de Sturm," Borel monographs, Paris, 1917. Rich indications 
as to the literature are found in the footnotes. 

f American Journal op Mathematics, Vol. XXXVIII, No. 1 (1916), particularly (1), (2), I, II, 
IV, V, VI and V. Sturm (Lioiiville's Journal,Vo\. I, p. 433), and Liouville (same Journal, Vol. I, p. 269) 
studied the interpolation problem, and the method of the latter has suggested that used here. But the 
only existing results appear to be confined to boundary conditions of Class I below, which does not include 
the periodic and other interesting cases. 

J Orthogonal here in the sense that C <j>i(x) <pj (x)r (m) das = 0, i^zj. 
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vanishes for x = for every i and j, and in consequence, by (3), also for x = l. 
We are thus led to the first class of boundary conditions, in which is also 
included the cases in which p(x) vanishes at one or both end points. 

Class I. M (0)^M(l)^(k One boundary condition involves only one 
end point, and the other, only the other. Consider first the case in which p (x) 
vanishes at an end point. Then under our assumption, all harmonics are 
different from at that end point. Consider an end point at which p(x) does 
not vanish, say x=0. Then, as M(0)=0, #(0)ft(0)— fc(0)#(0) =0, and 
since no harmonic vanishes simultaneously with its derivative, it follows that 
if any harmonic vanishes for x=0, all do. We conclude for Class I: If any 
harmonic vanishes at an end point, all vanish there. 

Supposing the determinant cd' — c'd = leads to the same class. We may 
therefore assume that both this and ab'—a'b are different from zero, and that 
p(0) >0 and p(l) >0. The conditions (2) may now be given the form 

y(l)=ay(0)+by'(0), (4 X ) 

y'(l) = cy(0)+dy'(0), (4 2 ) 

or y(0)= dy(l)-by'(l), (5,) 

y'(0)=-cy(l)+ay'(l), (5 2 ) 

where, by (3) and (4), A=ad — bc=p(0)/p(l) >0. These conditions charac- 
terize the second class of cases. 

Class II. M(0)=M(l)=fcO. The boundary conditions both involve both 
end points. M(0) and M(l) may both vanish, but not as a consequence of the 
boundary conditions. There will be three types to consider in this class. 

(a) b — 0. Then ad>0, and (4 2 ) shows that if $<(#) vanishes at either 
end point, it does at the other. The sign of d now becomes important. 

(a x ) &=0, a>0, d>0. If $i(x) vanishes at the end points, (4 2 ) shows 
that $4(0) and $(1) have the same signs, so that $i(x) must have an odd 
number of interior zeros, But (4^ shows the converse. Hence, for this case: 
All odd harmonics vanish, and all epen harmonics are different from zero at 
both end points. An example of this ease is y"-\-Xy — 0, y(l)=y(0), y'(l) 
=y'(0), with the .harmonics 1, sin 2nx, cos 2nx, sin inx, cos A-nx, 

(a 2 ) b — 0, a<0, d<0. By similar reasoning,. we conclude for this case: 
All even harmonics vanish, and all odd harmonics are different from zero at 
both end points. An example is y" + Xy=0, y(l) — — y(0), y'(l)= — y'(0), 
with the harmonies sin nx, cos itx, sin 3nx, cos Znx, 

(b) &>0. Since a harmonic will not vanish simultaneously with its 
derivative, (4 X ) shows that no harmonic vanishes at both end points. If 
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$>,.(0)=0, (4 X ) shows that <^(1) and^(O) have the same signs, and the har- 
monic must be even. The same conclusion follows from (5 2 ) when ^(1) =0. 
Hence in this case: No odd harmonic vanishes at either end point, and no har- 
monic vanishes at both. An example is y"-\-X 2 y=zO, ny(l)=2y'(0), 2y' (1) 

2 
= — 7ty(0). The harmonics are, for %>2, &(%) = siu \{x — 1) -| cos \x, 

where \ is the i-th positive root of 4^= (rf+AX 2 ) sin X; while \=X 1 —7t/2. 

For these frequencies we have the solutions A cos — x-\-B sin -x x. It will be 

seen if A and B are so chosen as to give an odd harmonic, this harmonic will 
be different from zero at both end points. All later harmonics are different 
from zero at both ends. 

(c) &<0. By similar reasoning we conclude : No even harmonic vanishes 
at either end point, and no harmonic does at both. 

3. Lemmas. 
We shall need the following : 

1. If two harmonics have the same frequency, the one with the greater 
number of interior zeros is different from zero at both end points. For the 
zeros of the two separate each other. 

2. Let y{x) =c q> Q (x) -}-c 1 ^) 1 (x) + . . . . -\-c n q> n {x) be a linear combination 
of harmonics, in which at least one of the constants c { before the last is not 0. 
If, for c n =0, y(x) has p interior roots, then for some c n ^0, it will have at 
least p interior roots. If for c n =0 it has an infinite number of interior roots, 
then for some c n =f=0, it will have at least N interior roots, N being any positive 
integer. The proof is easily supplied by considering the limit of y(x) as 
c„=0, the approach of y{x) to its limit being uniform. We shall reserve the 
notation y(x) for functions of the above form. 

3. Given a function y(x), c n =fz0, c n remaining fixed, we can choose 
c , c lt . . . ., c n _ x so small, not necessarily zero, that y(x) has the same number 
of interior zeros as $ n (x) (i. e., exactly n), provided <p n (x) does not vanish at 
an end point at which an earlier harmonic is different from zero. And in any 
case, y(x) has not more than one additional interior zero corresponding to a 
terminal zero of $ n (x). The proof depends upon the facts that $„(#) does not 
vanish simultaneously with its derivative, and that all the harmonics are finite, 
with finite derivatives. 
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4. Establishment of the Property (D). 

The results and treatment differ for the different types of boundary con- 
ditions enumerated. But the general argument is this. With a solution $>{x) 
of the differential equation, corresponding to boundary conditions to be 
assigned later, and not vanishing in the interior of (0, 1), we form the differ- 
ential operators to be applied to y(x) : 

Ny(x)=p(x)[y'(x)4,(x)—y(x)4,'(x)], and Ly{x) = ~Ny{x). (6) 

It will first be shown that repeated application of the operator L reduces 
relatively the earlier coefficients of y(x) to any desired degree, and thus, by 
the third lemma, leads to a function with no more than n interior roots. It 
will then be shown that the application of the operator L does not diminish 
the number of interior roots. The result, with the modifications noted, will be 
that no y(x) has more than n interior roots, and accordingly, that the deter- 
minant of the equations y(x ) =0, y(x x ) = 0, , y(x n ) =0 does not vanish 

for any set of arguments x { , no two of which are equal. The property (Z>) 
will follow (see this Journal, Vol. XL, p. 146, footnote). 

Note first, however, that D (x ) >0, since the first harmonic has no 
interior zero. Second, that in case ^i=A , Z>i(# , x y ) >0. For the two har- 
monics, belonging to the same frequency, will have a non- vanishing Wronskian, 
from which fact the desired result may be inferred. We may therefore sup- 
pose in the following, n>l, or, in case X X =X , n>2. 

The choice of -^(x) in the operators (6) may vary with the function y(x) 
to which it is to be applied, and the corresponding parameter value "k* for 
which ^{x) is a solution of the differential equation (1) will vary with it. We 
shall assume, however, that there are two constants, x and p, such that always 

x<l*<o where p< h+*2, or if ^=X , p<^-\ (7) 

The assumption will be considered below as use is made of it. Using the 
differential equation, it will be found that 

Ly(x)=p(x)^(x)r(x)[c (%* — 1 )<}> (x) +c 1 (%*— ^^(x) 

+ ....+c n (Z,*-\)<p n (x)]. 

Dividing, on the assumption c„=£0, by p(x)^{x)r{x) (A*— \) we arrive at a 
new function 

/^* ^ \ 
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/y* j^_ \ 

Under the assumption (7) the quantities (— s -1 are all less in absolute 

value than certain quantities independent of X* and less than 1, an exception 
arising only in the case of the last of them when X n = 2. _ 1 . Hence, by repeat- 
ing the application of L and the division, we may reduce all the coefficients of 
y(x) except the last (or, in case \—\_ x , the last two) to be less than any 
assigned positive quantity. This completes our first step. 

The next step is to study how this process affects the number of interior 
zeros in the various cases enumerated. 

Class I. Choose t^(x)=^ (x). Then 1*=7. , and (7) is satisfied. If 
y(x) has p interior zeros, Ny(x) will be seen to vanish between each pair, and 
thus have at least p — 1 interior zeros. But since in this case M(x) vanishes 
at both end points for all pairs of harmonics, and as y{x) is linear in the har- 
monics, it follows that Ny(x) also vanishes at both end points, and hence 
Ly{x) has at least p interior zeros. Thus no interior zeros are lost by the 
repeated application of L. Lemma 2 permits us to assume c n =f=0, and then 
Lemma 3 shows that no y(x) has more than n interior zeros, since by § 2, no 
harmonics vanish at an end point under the present boundary conditions 
unless all do. And a case of double frequency does not arise.* Hence, in 
problems of Class 1, the property (D) holds. 

Class II, (a^. & = 0, a>0, and y(l)=ay(0), for y{x) is a linear homo- 
geneous combination of harmonics, and hence satisfies the same boundary 
conditions. Hence y(0) and y(l), if not both zero, have the same signs. We 
again use -^(x) =q> (x), and assume that the last harmonic in y(x), $ n (x) is 
even. If y(x) has an odd number, 2p — 1, of interior zeros, it also has two at 
the end points, and Ny(x) has therefore at least 2p interior zeros, and Ly(x) 
at least 2p — 1. Thus, if the application of L to y(x) changes this function 
from one with an odd number of interior zeros to a new one with an odd 
number, none are lost. If it changes y (x) to a function with an even number 
of interior zeros, at least one is gained. If y(x) has an even number, 2p, of 
interior zeros, Ny(x) has at least 2p — 1, but (3) shows that it then either has 
an additional interior zero, or else two at the end points, and in either case 
Ly(x) has at least 2p — 1 interior zeros. If Ly(x) has an even number of 
interior zeros, it must be at least 2p, and if odd, at least 2p — 1. But ulti- 
mately, by Lemma 3, the repeated application of L leads to a function with an 
even number of interior zeros, since we have supposed this to be the case with 

* For two harmonics of the same frequency, M(x)z= const. This constant is not zero if the har- 
monics are independent. But it must be in the problems of Class I. 
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ty n (x). Hence no zeros have been lost, and we conclude: The harmonics of 
Case II {a x ) have the property (D) for every even n. 

The conclusion is still justified in case /l n _i,=/l n . Here we must study 
c n ^ n - 1 (x)+c n p n (x). If this has an even number of interior zeros, their 
number can not exceed n, since they are separated by the zeros of ^ n _ 1 (x) , and 
c„_i4>»_i(a0 +£„$>„(#) will, by the boundary conditions, be different from at 
the end points. The conclusion thus subsists in this case. If c n _ x q> n _ x {x) 
-\-c n q> n {x) has an odd number of interior zeros, the repeated application of L 
to y{x) may have reduced the number of interior zeros of y(x) by one, but by 
no more, and c^^^x) -\-c,$ n {x) has n — 1 interior zeros, and none at the end 
points, since by § 2, q> n -i(x) vanishes at both end points, and q> n {x) at neither, 
and e n 4=0. The application of the lemma 3 then shows that y(x) can not 
have more than n interior zeros, and the conclusion is established also in the 
case % n _!=\. 

Class II, (a 2 ). &=0, a<0, and y{l) =ay(0), so that if y(l) and y(0) 
are different from zero, they have opposite signs. We now operate with L in 
which ^(x) is chosen subject to the condition that y'{x)^{x) — y(x)4 , '{x) 
vanishes at both end points. By § 2, <p { x )> beiug an even harmonic, vanishes 
at both end points, and as ^(x) is different from zero on the interior of (0, 1), 
it follows that Pi* </l . Hence the assumption (7) is justified as far as the 
upper limit on Pi* is concerned. The existence of the lower limit will be 
generally established at the end of the paper. We start with a function y(x) 
ending in an odd harmonic. The operator N does not reduce the number of 
interior zeros by more than one, while with the present choice of ^{x), Ny(x) 
has zeros at the end points. So the operator L does not diminish the number 
of interior roots. As the odd harmonic <p n (x), by §2, is different from at 
the end points, Lemma 3 gives the result: The harmonics of Case II (a 2 ) have 
the property (D) for every odd n. 

Cases II, (b) and (c). Here, again, we choose ^(x) so as to make Ny(x) 
vanish at the end points. It is not, in these cases, apparent that the upper 
limit for PI* in (7) holds. However, since ^(^O does not vanish within (0, 1), 
and as q> % (%) vanishes twice, it is clear that Pl*<Pl a . As the frequencies have 
no finite limit point, there must be a positive integer n, and a number p, inde- 

pendent of Pi*, such that Pi* < p < °~T " . The least value of n for which these 
inequalities are possible,* we will call «'. The type of argument used in the 

* This least value of n can often be easily determined in a particular problem, it being helpful to 
observe that the value of X for which the differential equation has a solution vanishing at the end points 
of the interval (0, 1) but not in the interior, is an upper bound for X*. In the example in §9 under Case 
II (b), (7) is fulfilled without alteration. 
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previous cases then permits us to conclude: The harmonics of Case II (b) 
have the property (D) for every odd n>n', and, provided none of the har- 
monics after <£>„,_! (x) vanishes at either end point, for every n > n'. Also, the 
harmonics of Case II (c) have the property (D) for every even n>n', and, 
provided none of the harmonics after 4v_i(#) vanishes at either end point, 
for every n > «'. 

5. Existence of a Lower Bound for %*. 

In the problems of Class II, from (a 2 ) inclusive, on, it was assumed that 
a lower bound for the values of a,* existed. We proceed to show that such is 
the case. It will be recalled that ^{x) was to be chosen so that the ratios 
t^' ( 1 ) /^ ( 1 ) and ^'(0)/'4'(0) coincided with the corresponding ratios for y(x). 
Let us call these ratios v\ and £ respectively. They are not independent, but 
are connected by the equation, obtained by dividing (i 2 ) by (4 X ) : 

V=(c + d£)/(a + b£), where ad— bc>0. (8) 

We shall now prove the theorem : 

There exists a number A, independent of £ and y\, such that the solution of 
the differential equation (1) which satisfies the boundary conditions 

u'(0)=£u(0), u'{l)=m{l), (9) 

subject to the relation (8), and which does not vanish in the interior of (0, 1), 
belongs to a parameter value /1 # >A. 

We assume that p(0) >0 and p (1) >0. The theorem is not needed in the 
problems of Class I. 

We first prove the lemma : For fixed £ and decreasing yj, and for fixed y\ 
and increasing £, the value of the parameter a, corresponding to (9) increases. 
Let u, x (x) and u 2 {x) be two solutions of the differential equation satisfying the 
first equation (9). Denoting the corresponding parameter values by l x and \ , 
we find the following identity : 

p(l)[u 2 (l)u[(l)— u 1 (l)u' a (l)] = (% t — ^i) f u 2 (x)u 1 (x)r(x)dx. 

We are, in this paragraph, concerned only with solutions of (1) which do not 
vanish in the interior of (0, 1), so that it is legitimate to consider Wi(cc) and 
u 2 (x) positive for 0<«<1, so that the integral is positive. Dividing by 
u 2 (l)u 1 (l), we have 

p(l)[ yi (^)-y l (^))^ J^=^L C\(x)u 1 (x)r(x)dx, (10) 

U 2 (X)U 1 {1)^0 
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which shows that decreasing y\ corresponds to increasing A. The second part 
of the lemma is proved similarly. 

We next find the relation between £ and q when these are the values of 
the ratios u'(0)/u(0) and u'(l)/u(l) corresponding to any fixed "k. To do 
this we take a fundamental solution set for x=0, u x {x) and u 2 (x), for which 
1^(0) =1, 1^(0) =0, u 2 (0) =0, ^(0) =1, in terms of which any solution u(x) 
for the same value of a, may be expressed in the form u{x) —c x u x {x) +c 2 v 2 (x). 
Then % = c 2 /ci, and *i=[c 1 u' 1 (l)+c i u' 2 (l)]/[c 1 u 1 (l)+c !i u 2 (l)']. Eliminating 
the constants, we find the relation 

_ <(i)+gi4(i) nn 

^~«i(i)+5tf,(i)' K ' 

for which we have the easily derived relation 

u[(l)u 2 (l)-u 1 (l)u' 2 (l)=p(0)/p(l)>0. 

This shows that the hyperbola which is the Cartesian representation of the 
equation (11) is an always rising one. For various values of \ this hyperbola 
varies in position and size, but if for a finite value A of %, its center lies above 
and to the left of the hyperbola (8), also an always rising hyperbola, the 
lemma shows that the solution of the differential equation and the boundary 
conditions (9) subject to (8) will always correspond to a value of /1>A, and 
the theorem will be proved. 

That for some finite X the center of the hyperbola (11) lies above and to 
the left of the upper left-hand branch of the hyperbola (8), will be evident 
from the fact which we shall next prove : given any finite positive quantities 
h and k, a finite A, always exists for which the coordinates of the center of the 
hyperbola (11) satisfy the inequalities 

u 1 (l)/u 2 (l)>h, u' 2 (l)/u 2 (l)>k. (12) 

u (i) r 1 p(0) 

The identity ■ , ' = I ■ . . /. . dx, the derivation of which gives no 
1^(1) Jo p{x)u\{x) 

trouble if the initial values of the fundamental solution set are kept in mind, 

shows that the ratio %(!)/%(!) can be made as great as we wish if u^x) can 

be made as great as we please at all interior points. But this fact can be 

inferred from the identity 



d lu x {x, X)\ _ 



, ?V*' , r X f u 1 (x,X)u 1 (x,k 1 )r(x)dx. 



dx \Uxix, %x)/ p{ 

This shows that u x (x, ^)/m 1 («, \) increases with x if k<\. This ratio 
30 
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approaches 1 as #=^0. Hence always, u x (x, %) >u 1 (x, X) for 0<#<1, and 

using the law of the mean, 4- Ul( ; x >*\ > _kz^L C r {x)dx> {fa—X) X a posi- 

dx u-l{x, X) max. p J o 

tive function of x independent of X Hence the ratio, and therefore u x (x, X) 

can be made as large as we please, and the first inequality (12) is established. 

To prove the second, we start by noticing that for sufficiently small X 

u' 2 (x)>0 for 0<«<1. This is seen from the equation p(x)u' 2 (x)=p(0) 

+ f {—hr—q)u 2 dx, in which the integral is positive if /l<— max. q/m'm. r. 
Jo 

Multiplying the differential equation for u 2 by 2p(x)u' 2 (x) and integrating, we 



have [pu' 2 (x)]'' 



= — 2 J p(%r+q)u 2 u' 2 dx—— p(Xr-\-q)u\{l), the bar denoting 

n x'n 



1 



a mean value. Hence u 2 (l)/u 2 (l) > — ryr V(— A,) min. pr— max. pg, which 

can be made as large as we please by sufficiently diminishing X Thus the 
second inequality (12) is proved. 

Columbia, Mo., October 20, 1917. 



